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Abstract 

We study the thermodynamical properties of Fermi vapors confined in a 
harmonic external potential. In the case of the ideal Fermi gas, we compare 
exact density profiles with their semiclassical approximation in the conditions 
of recent experiments. Then, we consider the phase-separation of a multi- 
component Fermi vapor. In particular, we analyze the phase-separation as a 
function of temperature, number of particles and scattering length. Finally, 
we discuss the effect of rotation on the stability and thermodynamics of the 
trapped vapors. 
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I. INTRODUCTION 



Recent experiments with dilute vapors of alkali-metal atoms in magnetic or magneto- 
optical traps at very low temperatures have renewed the theoretical study of Bose and Fermi 
gases. For atomic gases, the Bose-Einstein condensation has been achieved in 1995 [1-3] and 
the Fermi quantum degeneracy in 1999 [4,5]. 

The s-wave scattering between Fermions in the same hyperfine state is inhibited due the 
Pauli principle. It follows that at low temperature the dilute Fermi gas, in a fixed hyperfine 
state, is practically ideal. Nevertheless, the effect of interaction could be very effective for 
a Fermi vapor with two or more hyperfine states (components). In the recent experiment 
with dilute 40 K Fermions [4,5], to favor the evaporative cooling, a 40 K Fermi vapor in two 
hyperfine states (|9/2,9/2) and |9/2, 7/2)) is used. When the system is below the Fermi 
temperature, one hyperfine component is removed and it remains a trapped quasi-ideal 
degenerate Fermi gas. 

In this paper we analyze the thermodynamical properties of a many-component inter- 
acting Fermi vapor by using a harmonic potential that models the trap of recent experiment 
with 40 K [4,5]. These trapped Fermi gases are quite interesting because the quantum degen- 
eracy shows up not only in momentum space, as in uniform systems, but also in coordinate 
space. First, we investigate the accuracy of the semiclassical approximation for the study 
of the properties of the Fermi gas in the conditions of 40 K experiment. For the ideal Fermi 
gas we discuss the finite-temperature spatial and momentum distribution and derive, in the 
large impulse approximation, the dynamic structure factor. Then, we consider Fermi vapors 
with many components. Due to the nonzero interaction between different hyperfine states, 
one can obtain the phase-separation of Fermi components by varying temperature, number 
of particles and scattering length. Finally, we discuss the effect of an imposed rotation 
around a fixed axis on stability and thermodynamics of the trapped Fermi vapors. 
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II. IDEAL FERMI GAS 



For an ideal Fermi gas at thermal equilibrium, the average number of particles in the 
single-particle state \a>) with energy e a is given by 

Ua = e fte a -n) + i ' ^ 

where /i is the chemical potential and (3 = l/(kT) with k the Boltzmann constant and T the 
temperature. The average number N = J2 a n a of particles of the system fixes the chemical 
potential [6]. In the case of a harmonic potential U(r) = (m/2)(ufxf + u)\x\ + ^3X3) one 
finds the exact density profile by using the Eq. ([I]) and Hermite polynomials, from which 
one obtains the eigenfunctions nin2ri3 (r) of the harmonic oscillator 

n( T \ = V |0nm 2 n 3 (r)| 2 , , 

V ; ^ e /3^(wi(m+l/2)+o;2(n2+l/2)+W3(n3+l/2)-A t ) _L 1 ' V > 

rtiri2n3=0 

The exact density profile can be compared with the semiclassical one. The semiclassical 
approximation means that, instead of e a , one uses the classical single-particle phase-space 
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energy e(r, p) = E— + U(r). In this way one obtains the single-particle phase-space distri- 
bution 

W ( r 'P) = ^(e(r,p) 1 -M) + l - 

At zero temperature, i.e., in the limit (3 —>■ oo where ji = Ep (the Fermi energy), one 
has n(r, p) = O (Ep — e(r, p)), where is the Heaviside step function. At temperature 
well below the Fermi temperature Tp = Ep/k, the Fermions begin to fill the lowest avail- 
able single-particle states in accordance with the Pauli exclusion Principle (Fermi quantum 
degeneracy) [6]. 

Because of the Heisenberg principle, the quantum elementary volume of the single- 
particle phase-space is given by (2iTh) 3 , where % is the Planck constant. It follows that 
the average number N of particles can be written as 

N = J n{r, p) = J dh n(r) = J d 3 p n(p) , (4) 



where n(r) is the spatial distribution and n(p) is the momentum distribution. The finite 
temperature spatial distribution is given by 



with T the factorial function. In the limit of zero temperature, with /i = Ep the 
Fermi energy, the spatial distribution gives the Thomas-Fermi approximation n(r) = 
(2m) 3 / 2 /(67r 2 ^ 3 ) [Ep - f/(r)) 3/2 6 (E F - U{r)), where 9 is the Heaviside step function [7-8]. 
The Eq. (^) is the generalization of well-known formula for an ideal homogenous Fermi gas 
in a box that is exact in the thermodynamical limit [6]. Note that for \z\ < 1 one has 
f n (z) = J2iZi(~^Y +lzl /i n - Moreover, by using g n (z) = —f n (—z) instead of f n (z), one finds 
the spatial distribution of the ideal Bose gas in external potential. 

A comparison between exact and semiclassical results has been recently performed at 
zero [9] and finite temperature [10]. It has been shown that the semiclassical approximation 
is good for [kT ' /Uuj h ) » 1, where ujh = (^i^^) 1 ^ 3 , or, at a fixed temperature, for a large 
number N of trapped particles. Our calculations confirm this prediction. In fact, we find that 
when (kT/hu>H) > 1 there are no appreciable deviations between exact and semiclassical 
results. Instead, when (kT/hutn) < 1 some differences are observable, in particular for 
"magic" numbers of particles (N = 1,4, 10,20,35,56,84, ...) that correspond to a complete 
shell occupation of single-particle energy levels (see also [10]). The differences are reduced by 
increasing N showing that semiclassical approximation provides an excellent representation 
of Fermi distribution for a wide range of parameters. In Figure 1 we plot the exact density 
profile, obtained by a direct calculation of Eq. ([J), and that obtained with the semiclassical 
function (|5|) for (kT/hun) = 10 -3 and three "magic" numbers. In correspondence of the 
"magic numbers", the spatial density profile shows local maxima, which suggest a spatial 
shell structure. The magic numbers are particularly stable; in fact, for small variations of 
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(5) 



where A = (27r^ 2 /?/m) 1 / 2 is the thermal length and 




(6) 



the chemical potential fi the magic number N remains unchanged. The shell structure in 
the density profile is washed out by increasing the number of particles and is completely 
absent in the semiclassical approximation. 

As shown in the experiment [4] , for the quasi-ideal one-component Fermi gas of 40 K atoms 
confined in a cigar-shaped trap (u>i = u>2 = 860.80 Hz, u 3 = 122.52 Hz), the semiclassical 
formulas of an ideal Fermi gas in a harmonic trap agree with the available experimental 
data within the statistical uncertainty. The accuracy of the semiclassical predictions are not 
surprising because the experimental conditions are kT/{huJn) — 10 2 -10 3 and N ~ 10 6 . For 
the sake of completeness, in Table 1 we show some relevant quantities of the Fermi gas at 
different temperatures: the energy n(r), the 

density at the origin n(0) and the widths (p 2 ) 1 ^ 2 and (z 2 ) 1 ^ 2 of the Fermi cloud in the radial 
and axial directions, respectively. 

In addition to the spatial distribution n(r), also the momentum distribution n(p) gives 
an useful signature of the quantum degeneracy of the Fermi gas. In the semiclassical approx- 
imation, the finite temperature momentum distribution of an ideal Fermi gas in a harmonic 
trap is calculated from Eq. (Q) and reads 

At zero temperature, one recovers the formula found by Butts and Rokhsar [7], namely 
n(p) = l/(6n 2 h 3 )2 3 / 2 /(muj 2 H ) 3 / 2 (E F - p 2 /2mfl 2 Q(E F - p 2 /2m). Despite the spatial 
anisotropy of the trap, the momentum distribution of the Fermi gas is isotropic when the 
semiclassical approximation is used. Moreover, it is not difficult to show that, both in the 
exact and semiclassical cases, the momentum distribution n can be found from the spa- 
tial one n by the simple relation h (pi,P2,P3) = n (pi/(mu>i),p2/(muj2),P3/(muj 3 )) /(mw^) 3 , 
where p = {pi,P2,P3)- in general, the exact momentum distribution is not isotropic and it 
presents the same shell structures of the spatial one when N is close to a magic number. 
Finally, the Eq. (|7|) can be used to calculate the dynamic structure factor ^(q, E) of the 
Fermi gas. In the large impulse approximation [11] one finds the remarkable result 
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where q and E are the momentum and energy transferred by the probe to the sample. 
Recently, the dynamic structure factor of a trapped Bose-Einstein condensate has been 
measured using two-photon optical Bragg spectroscopy with a time of flight technique by 
which the number of optically excited atoms can be counted [12]. 

III. MULTI-COMPONENT FERMI VAPOR 

The problem of a dilute Fermi vapor with M hyperfine states (components) can be stud- 
ied by using the s-wave scattering approximation, the mean-field approach and semiclassical 
formulas. The spatial density profile rij(r) of the i-th component of a Fermi vapor can be 
written as 



ni{v) = ^3/3/2 I e V > )\ 



(9) 



where i = 1,2,...,M, //j is the chemical potential of the i-th component, and = 
4irh 2 aij/m, with the s-wave scattering length between i-th and j-th component (an = 0). 
Thus, the effect of the other M — 1 Fermi components on the i-th component is the ap- 
pearance of a mean-field effective potential. Note that in the limit of zero temperature, one 
finds the equations used by Amoruso et al [13]. The total energy of a Fermi vapor with M 
components is given by 

M M 

+ d " r U(r)m(r) + Yl I d " r 9iM r hj(T) , (10) 

i=l i<j=l 

where the first term is the kinetic energy Ek in , the second term is the exter- 
nal potential energy E ext and the third term is the interaction energy E int . At 
zero temperature the kinetic energy assumes the familiar Thomas-Fermi form Eki n = 
E J =i(6vr 2 ) 2 / 3 (3a 2 )/(5m)/rf 3 r[n l (r)] 5 / 3 . 
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We numerically solve the set of equations @ with a self-consistent iterative procedure. 
Recently, we have used a similar iterative procedure to study 7 Li Bose gas [14]. If the 
components of the Fermi vapor are non-interacting then they can occupy the same spatial 
region. Instead, if the interaction is strong enough (repulsive interaction) or for iV very large 
one finds a phase-separation, i.e., the Fermi components stay in different spatial regions. This 
effect, that has been first discussed at T = by Amoruso et al [13], is shown in Figure 2, 
where we plot the density profiles of two components of the K Fermi vapor for different 
values of the scattering length a (a=ai2) and iV = 0.5 • 10 7 atoms in each component. 

When two components have the same number of particles, the onset of phase-separation 
is also an example of spontaneous symmetry breaking. In fact, if the chemical potentials /ij 
of the two components are equal, Eq. (|9]) always admits a symmetric solution ni(r) = n2(r). 
However, for particle number N larger than a threshold N c the solution bifurcates and a pair 
of symmetry breaking solutions appears. Just beyond threshold the asymptotic solutions 
begin to differ from the symmetric one in a neighborhood of the origin r = 0, i.e. at the 
point of higher density. Therefore, an analytic formula for the critical chemical potential /i c 
can be obtained by standard bifurcation analysis of Eq. ([]) at the origin, which gives 

At T = 0, by the use of the first term of the large z expansion of j 3/2 (^) [6], one finds 
analytical expressions for the critical density n c (0) and the critical chemical potential fi c : 

/~x 7T 57T 2 . f CLH\ 2 

nJO) = , u c = hu H — J • 12 

V ; 48a 3 ' P 24 V a J K ' 

These remarkably simple formulas can be very useful to determine the onset of phase- 
separation in future experiments. Moreover, by knowing the critical chemical potential one 
numerically finds the number of particles via Eq. (|^). 

At finite temperature, we numerically solve the Eq. (0). In Figure 3 we plot the 
critical number of Fermions as a function of the temperature for different values of the 
scattering length. Figure 3 shows that, by using a realistic scattering length [4,5], one needs 
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about N = 10 12 atoms to get the spontaneous symmetry breaking. One sees that at higher 
temperatures the phase-separation appears with a larger number of particles or a larger 
scattering length. In the lower part of Figure 3 we plot the same critical line as a function of 
temperature for different values of the scattering length. The conclusion is that by increasing 
the interaction between the two components one can use lower number of particles to obtain 
the phase-separation. Such a behavior can be extracted form Eq. (|12|) in the case of zero 
temperature. Note that at finite temperature one can use the first two terms of the large z 
expansion of jziiiz) [6]. Under the conditions kT << %ujh and \i >> gn(0) one finds the 
following equation for the onset of phase-separation 



One important point is to understand what happens when the two components of the 
Fermi vapor have a different number of atoms. In such a case, it is still possible to observe 
a phase-separation of the two components and their position is related to the ratio Ni/N%, 
where iV\ and N 2 are the number of Fermions in the two hyperfine states. In Figure 4 we 
plot the density profiles of the two components with different scattering length and ratios 
N1/N2. The component with less atoms is pushed outward but this effect appears only 
when the scattering length exceeds the value that produces the onset of phase separation 
with JVi = N 2 . 

Phase-separation also appears in a Fermi vapor with three or more components. In Figure 
5 we plot the density profiles of the 40 K Fermi vapor with three components. In this case 
we numerically solve Eq. @ with = GS13 — °23 = The Figure shows that also for three 
components the spontaneous symmetry breaking and the phase-separation are controlled by 
scattering length, temperature and total number of particles. In particular, one finds that 
by increasing the scattering length at first one of the components separates from the others, 
which remain still mixed. Note that the separation begins at the center of the trap. As 
the scattering length is further increased, also these two components separate (this second 
phase-separation begins at the interface with the previously separated component) and one 




(13) 



S 



eventually sees complete phase-separation and the formation of 4 or 5 shells. 

The Eq. ([11]) can be extended to a M-component Fermi vapor with the same number of 
particles in each component. The critical density n c (0) does not depend on the number M 
of Fermi components and one gets the same result of Eq. (pT2|). Instead, the critical chemical 
potential reads /i c = (2M + 1)tt 2 %ujh{cih I 'a) 2 /24. 

Nowadays it is possible to confine mixtures of Fermions and Bosons. If Bosons occupy 
the external region of the trap then the internal Fermi gas can increase its density at the 
origin and so favor the onset of phase-separation. We have verified that such effect is possible 
when two conditions are satisfied: 1) the number of Bosons is close to that of Fermions; 2) 
the Boson-Boson interaction is close to the Boson-Fermion one with scattering length of at 
least 10 3 Bohr radii. Nevertheless the gain is limited: the density at the origin grows by 20% 
for iV = 10 7 particles. In conclusion, to obtain the phase-separation with several millions 
of Fermions, the only chance is to strongly enhance the scattering length, for example by 
using the Feshbach resonances. The Feshbach resonances, already seen in Bose condensates 
[15], should allow to vary the interaction strength between atoms. One could use them to 
go smoothly from a quasi-ideal to a strongly interacting Fermi system. 



The semiclassical approximation can be easily applied to a rotating gas. For an ideal 
gas, that is rotating with angular velocity Q around the X3-axis, the classical single-particle 
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energy is e(r, p) = Jp- + U(r) — Q(xiP2 — x 2 pi). One immediately finds that the spatial 
distribution of an ideal Fermi gas is given by 



where p = (x\ + x 2 ,) 1 ^ 2 . Note that the effective potential of the rotating (fi > 0) Fermi gas 
is U(r) — (m/2)Qp 2 . If U(r) ~ p s with s < 2 for p — > oo then the effective potential goes to 
— oo for p — > oo and the spatial integral of the density profile is divergent. It follows that 



IV. ROTATING FERMI VAPOR 




(14) 
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for s < 2 the system is unstable against rotation. The same result applies for ideal Bose 
and Boltzmann gases. Obviously, the Fermi gas can remain confined near the center of the 
trap for a very long time depending on the tunneling probability. 

In the case of the harmonic potential, s = 2 and the gas is stable against rotation until Q 
is smaller than the smaller of the harmonic frequencies ui\ and uo 2 . The effect of the rotation 
is simply a shift in two trap frequencies: U\ — > (uj\ — Vt 2 ) 1 ^ 2 and 0J2 — > [lo\ — Vt 2 ) 1 / 2 . In this 
way one gets the Fermi energy of the rotating ideal gas 

^(") = 6'/W^(l-y . (15) 

A similar formula has been found in [16] for the Bose-Einstein transition temperature of 
a Bose gas in axial-symmetric harmonic potential. Also the density profile of the rotating 
quantum gas is modified according to the frequency shift discussed above. 

The discussion can be extended to a Fermi vapor with two or more components. In 
particular, one can use Eq. (|9]) and all the other equations by simply shifting the trap 
frequencies. Such a shift determines an enhancement of the characteristic length of the 
harmonic trap: a# — > clh/ (1 — fi 2 /co»i) 1 ^ 12 (l — Vl 2 /uj 2 ) 1 ^ 12 . It means that the rotating gas 
is more dilute and consequently the interatomic interaction is less effective. As shown in Eq. 
(|l2l) and Eq. (|13D , the rotation increases the critical number of Fermions that are necessary 
to get phase-separation. 



V. CONCLUSIONS 

First of all we have analyzed the accuracy of the semiclassical approximation. We have 
found that the semiclassical approach is very accurate for kT > Tiloh or, at a fixed tem- 
perature, for a large number of particles. In particular, the semiclassical formulas are ex- 
tremely reliable to study recent experiments with a trapped quasi-ideal 40 K Fermi gas where 
kT/(hujH) — 10 2 -10 3 and iV ~ 10 6 . By using the semiclassical approach, we have analyzed 
the spatial and momentum distributions of an ideal Fermi gas in harmonic potential. In par- 
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ticular, we have derived, in the large impulse approximation, the dynamic structure factor 
of the system. 

Then, we have studied the finite-temperature density profiles of a Fermi vapor with 
many hyperfine states. In the two-component case, by using the actual scattering length, 
we have found that the two components occupy the same spatial region. We have shown 
that the onset of phase-separation appears by increasing the scattering length or, for a fixed 
scattering length, by increasing the number of particles. By raising the temperature, a larger 
scattering length or a larger number of particles is needed to obtain the phase-separation. 
A Fermi vapor with three or more components has the same behavior but at first only one 
of the components separates from the others, which remain still mixed. The critical density 
of Fermions at the origin, which gives rise to the phase-separation, does not depend on the 
number of Fermi component and on the properties of the trap; moreover it satisfies the 
equation n c (0) = 7r/(48a 3 ), where a is the s-wave scattering length. 

Finally, we have considered the Fermi vapor is under rotation around a fixed axis with 
frequency Q. We have shown that, in the case of harmonic external potential, the rotation 
decreases the effective trap frequencies. As a consequence, the threshold for the onset of 
phase-separation occurs at larger numbers of Fermions. 
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T 


E/(kN) (p 2 ) 1 / 2 


(P 2 ) 1/2 


n(0) 


0.05 


0.44 


9.10 


44.60 


14.12 


0.27 


0.88 


12.89 


63.43 


7.37 


0.55 


2.66 


17.67 


87.26 


3.15 


0.82 


2.48 


21.56 


106.55 


1.76 


1.10 


3.30 


24.84 


122.90 


1.15 



TABLE 1. Some properties of the ideal Fermi gas with N = 7 ■ 10 5 atoms in the harmonic trap 
{u\ = LJ2 = 860.80 Hz, ix>3 = 122.52 Hz). Temperature T and energy per particle E/(kN) in units 
[xK, lengths in units fim and density n(0) in units 10 12 cm" 3 . 
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FIGURE CAPTIONS 



Figure 1: Comparison between exact (full) and semiclassical (dashed) density profiles 
for an ideal Fermi gas in an isotropic harmonic trap. kT/ihiOn) = 10~ 3 . iV is the number 
of Fermions. Lengths in units an = (h/muH) 1 ^ 2 and densities in units ajj 3 . 

Figure 2: Density profiles of the 40 K vapor with two components in the anisotropic 
harmonic trap. N = 0.5 ■ 10 7 atoms for each component, a) T = 0.1T F ; b) T = T F ; c) 
T = 2T F (Tp = 1.07 fiK). Units as in Figure 1. 

Figure 3: Critical number N of particles for the phase-separation of the 2-component 
40 K vapor as a function of the temperature T. On the top: a = 157a . On the bottom: 
a = 10 • 157ao (full diamond), a = 50 • 157a (open triangle), a = 100 • 157ao (full circle). a 
is the Bohr radius. 

Figure 4: Density profiles of the 40 K vapor with two components in the anisotropic 
harmonic trap at zero temperature. Ni = 10 7 atoms, a) a = 1000a ; b) a — 2700a (a is 
Bohr radius). Units as in Figure 1. 

Figure 5: Density profiles of the 40 K vapor with three components (solid, dotted and 
dashed lines) in the isotropic harmonic trap {uh = 450 Hz). iV = 0.5 ■ 10 7 atoms for each 
component, (a) T = 0.5T F ; (b) T = 2T F {T F = 1.07 /iK). Units as in Figure 1. 
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